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I. INTRODUCTION 


In many fields of physics, particularly quantum physics, atomic physics, subatomic physics 
and nuclear physics, for non relativistic case, no other fundamental equation has been studied 
more profoundly than the Schrodinger equation. Recently, considerable efforts have been made 
for several potentials to obtain analytical solutions of this famous equation, either in the ordinary 
case where the mass is considered as constant or in the context of position dependent effective 
mass. 

The position dependent effective mass formalism has been originally introduced by Von Roos in 
semiconductor theory ^ Later on, this formalism has been widely used in different fields of physics 
such as quantum liquids^, ^He clusters^, quantum wells, wires and dots^’^, metal clusters^, graded 
alloys and semiconductor heterostructures^"^^, the dependence of energy gap on magnetic field in 
semiconductor nano-scale quantum rings^"^, the solid state problems with the Dirac equation^^ and 
others^^“^^. Recently, it has been applied to study nuclear collective states within Bohr Hamilto¬ 
nian with Davidson potential and Kratzer potential^^”^"^. The advantage of this formalism resides 
in its ability to enhance the numerical calculation precision of physical observables, particularly 
the energy spectrum. Various methods have been used in the frame of this formalism for solving 
the Schrodinger equation for some potentials like the Point canonical transformation (PCT)^^“^^, 
Lie algebraic methods^^“^\ supersymmetric quantum mechanical (SUSYQM) and shape invari¬ 
ance (SI) techniques^^-^"^ and Nikiforov-Uvanov method^^. 

In the present work, following the precedure described in Ref. 20, we solve the Schrodinger equa¬ 
tion for some non central potentials by means of asymptotic iteration method (AIM)^^. For this 
purpose, we introduce a new generalized decomposition of the effective potential which allows 
simplifying the calculations through an easy separation of the coordinates. 

The AIM, an increasingly popular method, has proved to be a powerful, efficient and easily han¬ 
dling method in the treatment of problems in physics involving Schrodinger, Klein-Gordon and 
Dirac equations^^^. 

The content of this study is arranged as follows. In section II, we give basic concepts of the 
asymptotic iteration method. In section III, we present a theoretical background of the position 
dependent effective mass formalism. In section IV, the separation of variables is carried out for 
the Deformed Schrodinger equation (DSE) with general non central potential in spherical coor¬ 
dinates. In section V, we investigate the exact bound states solutions of Deformed Schrodinger 
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equation with different non central potentials. Finally, section VI is devoted to the discussion and 
conclusion. 


II. BASIC CONCEPTS OF THE ASYMPTOTIC ITERATION METHOD 


In this section, we present basic concepts of the AIM; for more details we refer the reader to 
Refs. 36 and 37. The AIM has been proposed to solve homogenous linear second-order differential 
equations of the form 


d^yn{x) , , .dyn{x) w W A 

= ?io{x)— - ^{x)tO 


( 2 . 1 ) 


dx^ ' dx 

Essentially, the functions 5o(.r) and Ao(.r) are sufficiently differentiable. Equation (2.1) can be 
iterated up to (fc-t-l)^^and (fc-f 2)'^ derivatives, = 1,2,3,---. Therefore we have 

d’^^^ynix) dyn{x) J*+2^„(x) dyn{x) ^ \ \ 

. k+\ =^k-\[x)— —+5fc_i(v)y„(;c), , ^+2 =h[x)— —+5,t(^)j/7(^) (2.2) 

ClJv CtJv CtJ^ 

where Xk{x) and Sk{x) are given by the recurrence relations 


Xk{x) = +2,o(v)Afc_i(v), -b5o(v)Afc_i(v) 


dx dx 

From the ratio of the {k-\-2y^ and {k-\- \ derivatives, we have 

d'^+^yn{x) Xkix) (-L 


(2.3) 


d , , (A:+l), 

LTiXyyi j - , J 

dx 




dx kk(x) 


yn{x) 


(2.4) 


-1^ + 

Now, we introduce the asymptotic aspect of the method. If we have, for sufficiently large k Refs. 
36 and 37, 


Skjx) ^ ^;^-l(jc) 
Xkix) Xk-l{x) 


:= a{x) 


(2.5) 


with a quantization condition 

. , , h{x) Sk{x) 

Afe(x) = 

Xk-\{x) Sk-\{x) 

then, the solution of Eq.(2. 1) can be written as 


= 0 , k= 1,2,3,--- 


(2.6) 


yn{x)=exp(-j a{z)dz 


C 2 + C 1 / exp 


(Ao(t) +2a{t))dt ] dz 


where Ci and C 2 are two constants. 

For a given potential, the procedure consists first to convert the Schrodinger equation into the form 
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of equation (2.1). Then, 5'o('*f) and Ao(a:) are determined, while Sn{x) and A„(a:) are ealeulated 
via the reeurrenee relations given by equation (2.3). The energy eigenvalues are then obtained 
by imposing the eondition shown in Eq. (2.6) if the problem is exaetly solvable. If not, for a 
speeifie prineipal quantum number n, we ehoose a suitable xq point, generally determined as the 
maximum value of the asymptotie wave funetion or the minimum value of the potential and the 
approximate energy eigenvalues are determined from the roots of the quantization eondition (2.6) 
for large values of k. As to the exaet eigenfunetions, they ean be derived from the following wave 
funetion generator: 

ynix) = C 2 exp J , k=\,2,---,n (2.7) 

III. FORMALISM OF POSITION-DEPENDENT EFFECTIVE MASS 


The aim of this seetion is to present a theoretieal baekground of the position-dependent effeetive 
mass formalism (PDEMF)^®’^^“^'^. In the PDEME, for Sehrodinger equation, the mass operator 
m{x) and momentum operator p — —ihV no longer eommute. Due to this reason and in order to 
obtain a Hermitian operator there are several ways for generalizing the usual form of the kinetic 
energy operator p^/2mo valid for a constant mass mo, since the generalization of the Hamiltonian 
describing the quantum state of a physical system is not trivial in this case. In order to avert any 
specific choices, one can use the general form of the Hamiltonian originally proposed by Von 
Roos^ : 


^ 4 L' 


(x)Vm'^(x)Vm^ {x)+m'^ (x)Vm'^ (x)Vm‘^ [x) -\-V{x) 


.k/ 


A'/ 


A' 


S', 


(3.1) 


where V (x) is the relevant potential and the parameters S', X', and k' are constrained by the con¬ 
dition : 5' + k' + X' ~ 

By choosing the position-dependent mass m(x) in the following form : 


m{x) = moM(x), M{x) = 


1 


/(^) 


(3.2) 


where mq is a constant mass, M{x) is a dimensionless position-dependent mass, and f{x) is a 
deforming function, the Hamiltonian in Eq. (3.1) transforms into 


// = — 


4 I 


r (x)Vm'^(x)V/^(x) +/^(x)V/'^(x)Vr (x) +y(x) 


(3.3) 


with d + K + X =2. The limit of the choice of the ambiguity parameters 5, K and X depends on the 
physical system. For example in the special choice where 5 = A = 0 and k = 2, the Hamiltonian 
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defined in Eq. (3.3) reduces to the most common BenDaniel-Duke form^'^. 

Following the method described in Ref. 20, the Hamiltonian in Eq. (3.3) can be written as 


// = -- 




2mo 


^/mynx)v^/m+Veff{x) 


with 




h 


2 r 


2mo 


i(1 -5-+ ('i-sVi-A) [Vf(x)f 




Therefore, the corresponding Deformed Schrodinger equation reads 


2mQ 


v^v/(x)vv7W+K//W 


\j/{x) = E\j/{x) 


(3.4) 


(3.5) 


(3.6) 


where E is the energy, h is the reduced Planck constant, t//(.r) is the total wave function and Veff{x) 
is the effective potential. 

IV. SEPARATION OF VARIABLES FOR THE DSE WITH GENERAL NON-CENTRAL 
POTENTIAL 


The problem of position dependent mass Schrodinger equation with a central potential has been 
widely investigated^^47,2i,35 unlike that with a non central potential. In spherical coordinates, the 
Deformed Schrodinger equation (3.6) for a particle in general non-central potential V{r,0, (p) reads 


as 




2mo 


with 


V f{r, 0, 9)V/(r, e,(p)Vxyf{r, e,(p)+Veff{r, 0, ^) ) V/(r,0,9) =Ev/(r,0,^) (4.1) 


(4.2) 


V,ff{r,e,<p) =V(r,e,(p) + ^|i(l-5-A)/(r,e,?>)vV(r,e,(p) 

In the special case of a mass {m{x) 1 /f{x)^) depending only on the radial variable r and in order 

to exactly separate the variables in Eq. (4.1), we propose a new generalized decomposition of the 
effective potential given by the following expression : 

f 1 


e, (P) = F(r, e, ?>) + ^ I j (1 - 5 - X)f(r)^-f(r) 

+ (i-A)lV/Wf 


(4.3) 
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where 


y(r,0,(p) = yi(r) + ^y2(0) + - ^ 


-y3((p) 


(4.4) 


' ' r^sirP-{0) 

with yi(r), y 2 ( 0 ) and y 3 ( 9 ) are arbitrary funetions of eertain arguments. 

This decomposition allows to introduce any forms of non-central potential in the Deformed 
Schrodinger equation (4.1). Besides, the particle wave function can be chosen in the follow¬ 
ing form 

1 R(r} 

(4.5) 




where the angular part of this function is selected in the form (p) = &{6)^{(p). 

Substituting Eq. (4.3) and Eq. (4.5) into Eq. (4.1) and using the standard procedure of separation 
of variables, one obtains the following equations: 


(fi ^ 2mo f E — Vi (r) 


dr^ 


+ 


/{rf 

1 


(2-5-1) ff"{r) ^ f\r) 


-8 


--A - 
2/4 


n ff(r) 


fir) 


R(r) = 0 




h 


0 ( 0 ) =0 


2mo 


d(p^ 




V3((p)+A^ 


d>(9) = 0 


(4.6) 

(4.7) 

(4.8) 


where we have introduced the seperation constants and L?' — £(£ + 1), where £ is the orbital 
angular momentum quantum number. 

V. SOLUTIONS FOR SOME NON-CENTRAL POTENTIALS 
A. Poschl-Teller double-ring-shaped Coulomb potential 


The Poschl-Teller double-ring-shaped Coulomb potential (PTDRSC) is defined by the follow- 
.,45 . 


mg expression 




1 


a^D(D-\) ^ a^C(C-l) ^ 


r \sirP-{9) ' cos^{0) J ' r^sirP-{9) V sirP-(a(p) ' cos^{a(p) 
whereA,C,D> l,/3 >0,b>0,a = 1,2,3,--. WhenA = C = D= 1 and= 0, the PTDRSC 
potential reduces to the Coulomb potential, which is one of the most important model in classical 
and quantum physics. When A = C = D = 1, the PTDRSC potential reduces to Hartmann poten¬ 
tial. Also, when C = D = 1, the PTDRSC potential reduces to the double-ring-shaped Coulomb 
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potential. 

Now, by adapting this potential to the general form given by equation (4.4) we have 


yi(r) 

V2(e) 

V3{(P) 


r’ 

b A(A-1) 
sin^{Q) cos^{Q) ’ 
a^D{D-l) a^C(C-l) 
sin^{a^>) cos^{cx<p) 


(5.2) 

(5.3) 

(5.4) 


Let us turn to find the analytieal solutions of the Deformed Sehrodinger equation with the PTDRSC 
potential. For this purpose and in order to find exaet analytieal results for Eq. (4.1) we are going 
to eonsider for the deformation funetion the speeial form^®’^^ : 


/(r) = 1 + er, e > 0 


(5.5) 


Generally, the ehoiee of the deformation funetion /(r) depends on the shape of the potential. 


I. Solutions of the radial equation 


When eonsidering the Coulomb potential shown in Eq. (5.2) and the deformation funetion 
given by Eq. (5.5), the radial equation (4.6) takes the form : 


2mo (Er-/3) \ e{2-5-X) ((1-2g) (1-2A)- 1) 

fp- r(l+er)2 \ r(H-er) 4(1+er)^ 


R{r)=Q (5.6) 


To solve this differential equation by means of the asymptotie iteration method, we propose the 
following ansatz 

;?(r) = r“(l+er)''g(r) 


with 


where 


1 + \/1+4L2 ^ _ 1 ^ ^£2(2M-l)2 + 4(£a-T) 

^ — "Z 1 ^ — T ^ Z 

2 2 2 £ 


ct = ^^ + (5+A-2(1+l 2))£, r=^!^+Q(5 + A)-(2+L2 + 5A) )£2 

Eor this form of the wave funetion, the radial equation (5.6) reads 

d'^gnXr) .dgnfr) , , , , , 


(5.7) 


(5.8) 


(5.9) 


dr^ 


dr 


(5.10) 
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with 


AoW = 


2£r{u + v) +2u 




2£u{u “ 1 “ V — 1 ) “ t “ c 


(5.11) 


r(l +er) r(l +£r) 

According to the AIM procedure, the energy eigenvalues are then computed by means of the 


quantization condition (2.6) : 


VO = 


Vl = 


V2 = 


V3 = 


2£{u^ — u)+o 
2£u 
2£u^ + o 
2£{u + 1 )’ 

2c (m^ “1“ m a 1) “t“ c 
2c (m -\- 2) 

2c(m^ 2u A3) c 
2c(m -|- 3) 


(5.12) 


which can be generalized to 




2C + {Ur —\)u + 




+ 0 


n, = 0, 1,2,--- 


(5.13) 


2c(M + n^) 

Substituting u and v by their expressions given in Eq. (5.8), we finally derive the exact eingenvalues 
of the radial equation : 


Enr,e — 


/^“^(^(^+l)+2-5-A)cj mo ;j 2 g 2 

^ {nr + (. + \Y 


2r {tir+i+iy 


8mo 


h 


+ ;^(/3 + ^(^(^+l) + 5 + A)c) + 




2mo 

n^ = 0,1,2,3,-•• 


1-5-A+ --5 --A c2. 


(5.14) 


which, in the c —0 limit, leads to the usual result = —mo/3^/(2^^(n^+ ^+ 1)^), where 
/3 = Ze^. The corresponding eigenfunctions are the hypergeometrical functions, 

gnXr) = iFi {-nr,nr + 2v + 2u- l;2v; 1 +cr) (5.15) 

Consequently, the radial wave functions reads as 

Rnr/ir) =C„^r"(l +cr)'' 2 Fi (-n^,n^ + 2v + 2M- l;2v; 1 +Cr) (5.16) 


where C„^ is a normalization eonstant. Using the normalization eondition of the radial wave fune- 
tion and the following series representation of the hypergeometric fucntions'^^’'^^ : 


p 


Fq (fill, • • • ■ ■ 


oo 

Fq\z) = Y, 


n=0 


(ai)n---(ap)nZ" 
(fil )n ■ ■ ■ Yq)n 


(5.17) 
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we obtain the normalization eonstant 


{rir + 2v + 2u — 


1). 


r(2M + i)-<2, 


(“>v; 

rir 


- 1/2 


with 


(5.18) 




^ (-n,)^(l-2v-n,)^ ^ (-n.); (l-2v-n,)^. r(fc + 7 +3 - 2v - 2 m - 2n,) 

^Qk\{2 — 2v — 2u — 2nr)k ^Qjl{2 — 2v — 2u — 2nr)j T{j+ k + 4 — 2nr — 2v) 

(5.19) 


where the Poehhammer symbols (a); are defined by 


{a)i 


nci + i) (-l)^r(n, + l) 

r(a) r(n,-A:+l) 


(5.20) 


2, Solutions of the first angular equation 


We are now going to derive eigenvalues and eigenfunetions of the angular equation (4.7). Using 
the potential given by Eq. (5.3), one ean write the differential equation (4.7) as 


(fi / d 2 ds?' 

^_+co,(e)-+L -- 


2mo f b A (A — 1) 
\ sin^{Q) 


0 ( 0 ) =0 


(5.21) 


sin^{6) \sin^{6) cos^(G) 

In order to apply the AIM approaeh for the above equation, we introduee a new variable y = 
cos(9) e [—1,1]. So, we obtain the following differential equation 

d^&(y) 2y d&(y) (-LV + icy^-y) 


dy^ 


1 — y2 (iy 


■ + ■ 


(y(l-y2))2 


-0(y) = 0 


with 


,2 , .2 , 2mo ,,,, 2mo 


k: = L^ + A^ + -^(A(A-1)-Z;), 7=-^A(A-1) 
h h 


A reasonable physieal wave funetion is proposed as follows : 


(5.22) 


(5.23) 


0W=3'’'(1-A'’«W 


(5.24) 


with 


0 




8/noA(A—1) 
¥ 


2 


P 



The equation (5.22) reduees to 


d'^^kjy) 

dy^ 


>^o(y) 


d^k{y) 

dy 




(5.25) 


(5.26) 
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with 




2((T7+2p + l)/-?7) 

3^(1 -y^) 


, ^o(j) 


((?7+2p + l)(T7+2p)-L2) 
1 — 


Thanks to the quantization condition (2.6), we derive the energy eigenvalues as : 


(5.27) 


4 = T?+2p, 

£i = r7+2p + 2, (5.28) 

£2 = t?+2p+4, 


from which, we get the general form of the quantum number i 

ik — V 2p + 2k, k = 0,1,2, - ■■ 


(5.29) 


Substituting rj and p by their expressions given in Eq. (5.25), we finally obtain the full expression 
of £ 


+ \/A^+ + 2k, k = 0,1,2,--- 


(5.30) 


If we set t = y^, the differential equation (5.26) transforms into the well-known Gauss hyper¬ 
geometric differential equation Eq.(15.10.1) in Ref.46. So, the solution of Eq.(5.26) is given by 


^{y) — 2p\ k,k-t-7] -t-2p-|--;7] + 2’^^^ (5.31) 

Einally, the wave function solutions of equation (5.21) are obtained : 

0(0) = Ck{cos{e))^{sin{e)fp 2 Fi (^-k,k + r] +2p + + ^^cos^e)^ (5.32) 

where Q is the normalization constant of the angular wave function 0(0). This constant is calcu¬ 
lated from the normalization condition : 


rsin{e)\&{e)\^de =2 C |0(y)|2jy= l 

Jo Jo 

By using the following relation of the orthogonality of the Jacobi polynomials 

^2 


46,47. 


Jo 

we obtain 


^ .s-Yi 17 / , m 2 , r(7)^r(n-t-5-7-f 1) 

1-z) ^[ 2 Fi{-n,n + s;r,z)] dz^ , , . , , , — 

[s + 2n) r(n-b5')r(5-b7) 


Ck = 


\ 


(2k -f T] 4- 2p -|-2)r(k-|-7]4- 2p 5 ) T (k -f p 4- 2 ) 


k\ 


r{r] + ^Yr{k+2p + i) 


where p and p are given in Eq. (5.25). 


(5.33) 


(5.34) 


(5.35) 
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3. Solutions of the second angular equation 


After introducing the expression of the potential given by equation (5.4) into (4.8) we get 




2 lmQfa^D{D-\) a^C{C-r 
V sin^{a(p) ^ cos^(a(p) 


= 0 


Defining a new variable z = cos{a(p), the differential equation (5.36) transforms to 
J 2 $(z) z d^{z) {-e^z^ + CJz^-Td) 


with 


dz^ 


2 ^ 

£ =^, dJ = 
a? 


1—z2 dz 


+ ■ 


(Z(l-Z 2))2 


2mo(C + D-l)(C-D) 


f e , -& = 


d>(z) = 0 
2moC(C-l) 


We choose the following ansatz 


d>(z)=z^(l-z^)''z(z) 


with 


1 


1 , 8moC(C-l) 


1 

■'’=4 + 


1 + 


8mo£)(£)-l) 

fr 


When we insert Eq. (5.39) into equation (5.38) we get 

72, 


d Xq{z) _ , , .dXqiz) 


dz^ 


= Ao(z)- 


dz 


'rSQ{z)Xq{z) 


with 


a,(.) = PM+4 v+1)^-2m ^ ^ (2v + m) 


2-£2 


Z(l-z 2 ) ’ l-z 2 

By using the same procedure as in the previous case and after few iterations we get 

Co = ±(2v + jU), Cl = ±(2v + jU+2), £2 = ±(2v + jU+4), £3 = ±(2v + jU + 6 ), 

which induce the generalized relation of £ 

£^ = ±(/i+ 2 v+ 2 ^), ^ = 0 , 1 , 2 ,- 

This equation can be rewritten in the following form 


A.q — ib(x(/i+ 2 v+ 2 ^), q — 0 , 1 , 2 ,--- 


Substituting the obtained expression for /i and v Eq. (5.40) into Eq. (5.45) we get 

^ 8moC(C-'TT A ^ 8mo£>(P^ \ 

A = ±a I + 2^+1 , ^ = 0,1,2,••• 


(5.36) 


(5.37) 


(5.38) 


(5.39) 


(5.40) 


(5.41) 


(5.42) 


(5.43) 


(5.44) 


(5.45) 


(5.46) 
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Like equation (5.26), the eigenfunctions of equation (5.41) are the hypergeometrical functions 


X{z) — iFi ( —q,li +2v + ^;/i + 


Finally, we find the wave function of the angular equation (5.36) 


=Cq{cos{a^)))^{sin{a(p)f-'^2F\ ( -q,pL+2v+q-,ii +^^cos^^a^)) 


1 


where Cq is the normalization constant of the angular wave function given by 


(5.47) 


(5.48) 


Q = 


(H+lv + 2 q) r (/J + 2v + ^) r (v + ^ + 2 ) 

\| Z'?! r(M + yV(2v+«+i) 


where /i and v are given in Eq. (5.40). 

B. Double ring-shaped Kratzer potential 

In spherical coordinates, this potential is defined as^^ : 

1 


y(r,0,(p) = -2D, 


+ ■ 


b a 

+ 


, a > 0, h > 0 


(5.49) 


(5.50) 


r 2r^ J \sirF-{d) cos^{d) 

where Dg is the dissociation energy between two atoms in a solid and rg is the equilibrium inter- 
nuclear length. In this cas we choose 

.2X / ^ 




. V3(<p) = 0 (5.51) 


To find the energy spectrum and eigenfunctions of the Schrodinger equation with double ring- 
shaped Kratzer potential, we need just to replace L?- = £{£ + 1) by = L^ + expres¬ 

sion of u in Eq. (5.8) and to set /3 = 2Dgrg in the expression of v in Eq. (5.8) : 


u- 2 + 




_1 , \/e^(2«-l)2+4(e(7-T) 

'^ 2 '^ 2e 


(5.52) 


with 


a = + {5 + X -2{1 +L^)) £, t = ^^+(1{5 + X)-{2+L^ + 5X) ) (5.53) 
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By using the identity relating m to v Eq. (5.13), we finally deduee the energy speetrum for double 
ring-shaped Kratzer potential, 


Enr,e —' 


2Dere-^^{2e{i+l) + 2-5-X)£y mo 




2 h^ {jlr + U 
-.2 


2mo{nr-\rU)^ 


nr{nr+\){u + ^) 


n — 1 2Pi^f ^ \ ) 

X + + - -u{u-\) 2D,r,-;—(2£(£+l)-l-2-5-A) \ 

2 mo V 2mo ) J 


+ 


2mo 


{l-5){l-A,)£+Dere e, n^ = 0,1,2,--- 


(5.54) 


where the quantum number £ is dedueed from equation (5.30) : 


1 +Jl + ^ 


-f A/a^- h -f 2fc, fc = 0,l,2,--- 


(5.55) 


In addition, the radial and angular wave funetions are identieal to those given by the equations 
(5.16) and (5.32) respeetively. On the other hand, from the results obtained for Kratzer potential 
(KP) we ean also deduee the eigenvalues and eigenfunetions for a new type of potential ealled the 
Modified Kratzer potential (MKP) or Kratzer-Fues potential defined as"^^’^^ : 

2 


Vi{r)=De 


r-r. 


(5.56) 


whieh is obtained by adding a Dg term to the potential Vi (r) in Eq. (5.51). In this ease, the energy 
speetrum formula is simply, 

(5.57) 




"nr,£ 

(KP) 

where E « is the energy spreetrum given by the formula (5.54). 


C. Makarov Potential 


This potential has the following form^^’^^ : 


V(r,e,tp) 


r \sin^{6) sirfi{Q) J ’ 


(5.58) 


The first term is the Coulomb potential, the seeond and the third represent the short-range ring- 
shape terms. The Makarov potential ean be used to deseribe ring-shaped moleeules sueh as ben¬ 
zene and interaetions between deformed pairs of nuelei. The radial energy speetrum and the radial 
wave funetions eorresponding to this potential have been already determined in the subseetion 
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(A.l). So, we shall study only the analytieal solutions of Eq. (4.7). For this purpose, we eonsider 
the expression of the potential ¥ 2 ( 6 ) as 



sin^{9) J 


(5.59) 


Insering this potential into equation (4.7), we obtain the angular equation as 


d 

+ cot( 0 )— +L' 




2mo f Ct 


<702 ' ' sin^{6) \sin^{6) sirP-{9) j 


Ycos(O)) 


0 ( 0 ) =0 


(5.60) 


We define a new variable t = 5(1 +cos(6)), so Eq. (5.60) beeomes 
d^&(t) (l- 2 t)d&(t) (-pV + ut + to 


dt^ t(l — t) df 
where the eonstants are given by 


t 2 (i -ty 


0(0 = 0 


p=L, v=L^- ^7, (0 = ^ (r- a) - ^A^ 

2h^ ^ 4 


(5.61) 


(5.62) 


We take the wave funetion in the form 


0(0 = t“{\ -tyh{t) 


(5.63) 


where 


“ = + v = ^jA2 + ^(a + 7) 


(5.64) 


2 V r ' 2 V 

By inserting the ansatz given in Eq. (5.63) into Eq. (5.61), we obtain the following differential 
equation 


d^hj{t) 

(iO 


= Ao (0 


dhj{t) 


dt 


¥sQhj{t) 


with 


2(m + v+1)^-2m-1 (m + v+1)(m + v)-^2 

'^(O =-—A-’ ■^o(0 = — 


(5.65) 


(5.66) 


01-0 ’ ^( 1-0 
Aeeording to the AIM proeedure, the energy eigenvalues are then eomputed by means of the 
quantization eondition (2.6). After few iterations we obtain 


^Q = U-\rV, ^\ = U + V £2 = M + V + 2 , £3 = M + V + 3 , ••• ( 5 . 67 ) 


In general form, we have 


ij = u + v + j, j = 0,l,2,--- 


(5.68) 
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Substituting u and v by their expressions given in Eq. (5.64), we finally derive the exaet formula 
of the quantun number ^ 


= ^yA^ + ^(a-7) + ^YA2 + ^(a + 7) + y, ; =0,1,2, 


+ 7 , 7 = 0 , 1 , 2 , 



(5.69) 


The Eq. (5.65) is the hypergeometrie differential equation, whose solutions are given by 


h{t) = 2^1 (-7,2m + 2v + 7 + 1;2m + l;t) 


(5.70) 


Einally, the angular wave funetion (5.60) ean be written as 
1 +co5'(0)\ “ /1 —cos(O) 


0 ( 0 ) =c,- 


^ f . . l+cos(6) 

2^1 ( ~7; 2m + 2v + 7 + 1; 2 m + 1;--- 


(5.71) 

where Cj is a normalization eonstant computed via the orthogonality relation of Jacobi polynomi¬ 
als 


C; = 


'{u + v + j + j) r {2u + 2v + j + l)r {2u + j + 1) 


7! 


r{2u + irr{2v+j+i) 


(5.72) 


where u and v are given by the formula (5.64). 


D. A new Coulomb ring-shaped potential 


In this subsection, we solve the angular part of Eq. (4.7) with two potentials using AIM . 
The first potential Vie(6) is a novel angle-dependent (NAD) potential, originally introduced by 
Berkdemir^®’^^ and the second potential V2e{0) is another kind of NAD potential, introduced by 
Zhang and Huang-Eu^^. These potentials are defined as 

f r+Ksin^{e)+r]sin\e) 

V sin^{0)cos^{e) 


(5.73) 


V2e{e) 


Y+ kcos^{ 6) + ricos^{6) 
sin^{9)cos^{0) 


(5.74) 


where 7 , k and rj are real constants. 

The last potential is deduced from the potential in Eq. (5.73) using a simple transformation of 
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the angle 6 : ¥ 20 ( 6 ) = Vie{0 ^ 0 ± f). So, in our case we can express the NAD-type Coulomb 
potential as 


V{r,e,cp)={ 


jS I ft" ViejO) 

r ' 2mo 


I fP' ViejO) 
r ' 2m{) 


(5.75) 


where a = is the fine structure constant and /3 = Zct. 

The radial energy spectrum and the radial wave functions of the standard Coulomb potential are 
already similar to those obtained in the previous section. Let us now turn to the calculation of the 
energy spectrum and the normalized wave functions for the first potential. The substitution of the 
above potential (5.73) into Eq. (4.7) leads to the following differential equation 

^ Y-^Ksirp-{d)+r\sin^{d)' 


d'^ , d 9 

—+com—+L 


0 ( 0 ) =0 


sin^(6) sin^{0)cos^{6) 

To solve this equation, we first introduce the transformation x = cos^{0), so we get 


(5.76) 


d^&(x) ( 1 —3x) d&(x) (—px^ + vx—K) 

41-4 


(5.77) 


with 

p = ^{L^ + j]), v = ^{L^ + k + 2j]-A^), K = ^{Y+r] + K) (5.78) 

In order to obtain an equation suitable for applying the AIM approach, we choose the wave func¬ 
tion as 


&{x) = x‘i{i-x)Pi;{x) 


(5.79) 


with 

P = ^ + ^Vl+4(7+T] + ?f), q = ^^A2 + Y 
Inserting it into Eq. (5.77) we get 


d^Cijx) 

dx^ 




(5.80) 


(5.81) 


with 


.. {p + q + l){P + q)-H^^ + Pl) {2p + 2q+^)x-{2p + ^) 

=- 41 ^^^- 


(5.82) 
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Applying the AIM to Eq. (5.81), the identities that eonnects p and q are then eomputed by means 
of the quantization eondition (2.6). For few iterations, the obtained solutions are : 


— -(<?0 + ±^V^1+4(L2 + 7]), 

Pi = — (^1 + 1 +-) ± - (5.83) 

P 2 = -(^2 + 2 +^)±^ y / i + 4 ( L 2 + 77 ), 

Then we deduee the general formula : 

p — — ±^"y/l +4(L2 + 7]), z = 0,1,2,• • • (5.84) 

From the last relation, we obtain the formula of 

= 4+ (8 (p + ^) + 2)/ + (8p + 2)<5'+ 2p — 7], z = 0,1,2,--- (5.85) 


By insering the values of p and q given by Eq. (5.80) in Eq. (5.85) we finally get 


L 2 = l+4z(i+l) + ?c + 27 +A 2 + 2 ( 2 z + l)^A 2 + 7 

+ ^1 + 2z + \J A2-|- 7 ^ a^1+4(7+77 ~ 0,1,2, • • ■ (5.86) 

Essentially, the eigenfunetions of Eq. (5.81) are the hypergeometrieal funetions 

l^{x) =2^1 |^-z,z + 2p + 2<5r + ^;2p + ^;A:^ (5.87) 

From equations (5.87) and (5.79) we obtain the normali z ed wave funetion 

0(0) = Q {cos{d))^’' 2 F 1 (-i, i + 2p + 2q + ^;2p + ^;co 5 '^( 0 )'j (5.88) 


with 


Ci = 




(p + ^ + z + ^) E {2p + 2 ^ + z' + E {2p + * + 5 ) 


Z! 


E(2p + i) E(2^ + z + l) 


(5.89) 


where the parameters p and q are given in Eq. (5.80). 

On the other hand, if we insert the potential VzeiG) (5.74) into Eq. (4.7), then after some algebraic 
manipulations, we obtain the same differential equation (5.81) where the parameters p and q are 
taken as 

P = ^ + ^V1+47, ^V^A2 + 77 + k + 7 (5.90) 
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Hence, the obtained energy spectrum and wave functions in this subsection for the novel angle- 
dependent potential Vie{6) are also valid for the NAD potential V2e{0)- Consequently, Eq. (5.85) 
turns to 


= 1+4/(/+ 1 )-f-K-h 27 +a2 + ( 2 /+ 1 ) 



(5.91) 


which is consistent with the result reported in Ref. 53. 

VI. DISCUSSION AND CONCLUSION 

In this work, by means of AIM, we have solved the Schrodinger equation, within the for¬ 
malism of position dependent effective mass, for a class of non central physical potentials such 
as Pdschl-Teller double-ring-shaped Coulomb potential, Makarov Potential, double ring-shaped 
Kratzer potential, double ring-shaped Modified Kratzer potential (Kratzer-Fues potential) and a 
new Coulomb ring-shaped potential (Coulomb potantial plus novel angle-dependent potential). In 
order to simplify such calculations, we have introduced a new generalized decomposition of the 
used effective potential. In our approach we have shown that the used new decomposition of the 
effective potential allows studying DSE with different non-central potentials. Moreover, we dis¬ 
cuss some special cases of interest which can be deduced from our general solutions obtained in 
this work. So, we start with the first special case where £ 0 corresponding to DSE with a purely 

central potential: in the absence of the ring-shaped form interaction, the used non central potential 
reduces to standard Coulomb potential and the corresponding energy spectrum given in Eq. (5.14) 
is identical to the result obtained in Ref. 20 via supersymmetric quantum mechanical and shape 
invariance techniques. While, in the limit £ —> 0 corresponding to standard SE with non central 
potentials, our results, for the used type potentials, reproduce exactly those previously obtained 
in several works especially for double ring-shaped Kratzer potential'^^, Pdschl-Teller double-ring¬ 
shaped Coulomb potential'^^, Makarov Potential^^, a new Coulomb ring-shaped potential^^ and 
double ring-shaped Coulomb potential^^. The present work could be extended to generalize the 
solution in multi-dimensional space of Deformed Schrodinger equation with any non central po¬ 
tential. Finally, the obtained theoretical results could find many applications in several fields of 
physics. 
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